Astronomy & Astrophysics manuscript no. CO vib AuA 


© ESO 2009 


February 19, 2009 





o 
o 

(N 
0\ 



a 

Oh; 
O 



> 
00 

m 
m 

(N 
O 

o\ 
o 



Accurate laboratory rest frequencies of vibrationally excited CO up 

to i; = 3 and up to 2 THz 

R. Gendriesch^ F. Lewen\ G. Klapper\ K. M. Menten^, G. Winnewisser\ J. A. Coxon^, and H. S. P. Muller^ 

' I. Physikalisches Institut, Universitat zu Koln, Zulpicher Str. 77, 50937 Koln, Germany 

e-mail: lewen , winnewisser , hspm@phl . uni -koeln . de 
- Max-Planck Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

e-mail: kmentenOmpifr-bonn.mpg . de 
^ Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, B3H 4J3, Canada e-mail: John. Coxon@Dal . Ca 

Received 18 December 2008 / Accepted 06 February 2009 

ABSTRACT 

Aims. Astronomical observations of (sub)millimeter wavelength pure rotational emission lines of the second most abundant molecule 
in the Universe, CO, hold the promise of probing regions of high temperature and density in the innermost parts of circumstellar 
envelopes. 

Methods. The rotational spectrum of vibrationally excited CO up to o = 3 has been measured in the laboratory between 220 and 
1940 GHz with relative accuracies up to 5.2 x 10"', corresponding to ~ 5 kHz near 1 THz. 

Results. The rotational constant B and the quartic distortion parameter D have been determined with high accuracy and even the 
sextic distortion term H was determined quite well for v = I while reasonable estimates of H were obtained for v = 2 and 3. 
Conclusions. The present data set allows for the prediction of accurate rest frequencies of vibrationally excited CO well beyond 
2 THz. 
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1. Introduction 



methods: laboratory - techniques: spectroscopic - radio lines: ISM - radio lines: stars - ISM: molecules 



Carbon monoxide, CO, is the second most abundant molecule 
in the universe next to hydrogen, H2. Observations of CO are 
commonly employed to obtain information on the H2 abun- 
dance since H2 is usually difficult to detect directly because it 
does not have any permanent or transition dipole moment, and 
its quadrupole transitions are very weak and high in energies. 

^ Usually transitions of the main isotopologue '^C'^O (in the fol- 
lowing, unlabeled atoms designate '^C or '^O) are recorded for 
this purpose, but frequently also transitions of '^CO or even 
C'**0 are studied because transitions of these isotopic species 
are less affected by opacity. Even less abundant isotopologues 

' have been de tected in space, inc luding recenfly the least abun- 
dant i^c'^O jBensch et al.ll200ll) . 

Very accurate ground state rest frequencies have been ob- 
tained in the laboratory for CO jW innewisser et al. 1997) as we ll 
as for its isotopologues (IZinketal.„1990c.Klapi7er et al.l2000allbl 
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120031; ICazzoU et al.1 l2002l l2003l l2004t iPuzzarini et al.1 
20031) . The lower-7 transitions were usually measured in sub- 
Doppler resolution such that the relative accuracies of the mea- 
surements reach ~ 10"''. 

Interestingly, the only experimental article to report mea- 
surements of transitions of vibrationally excited CO seems to 
be a paper by iBogev et al.l (Il986l) on CO 7 = 2-1 transi- 
tions for several vibrational states between u - 5 and u - 
40. One reason for this lack of measurements is probably the 
high vibr ational energy of 2143 cm"' or 3086 K for u - I 
dCoxon & Hajigeorgiou 2004) which makes excited states of CO 
difficult to populate. Anot her reason is l ikely the small dipole 
moment of 0.10980 (3) D (iMuenteill975h which even decreases 



for the first few vibrational states by about 0.025 D for each ad- 
ditional vibrational excitation iintil the dipole m oment switches 
polarity between v - 4 and 5 (lGoorvitchlll994l) . see also Fig[I] 
ICoxon & Haiigeorgioul ( |2004|) reviewed rather extensively not 
only rotational data but also a plethora of infrared studies. Their 
work permits very good predictions of the rotational transition 
of CO in excited vibrational states. 

Apart from molecular clouds, CO is also a premier con- 
stituent of the circumstellar envelopes of mass-losing stars on 
the asymptotic giant branch (AGB) of both carbon-rich ([C/O] 
abundance ratio > 1) and oxygen-rich ([C/O] < 1) objects. 
Moreover, unlike other species, it does not condense into dust 
grains and prevails throughout the envelope until it is destroyed 
by ultraviolet radiation in its tenuous otiter po rtions (see, e.g., 
iLafont et al.lll982t ICherchne ff & Barkeii ll992h . Thus, emission 
from its strong rotational lines dominate the cooling and, thus, 
the th ermal balance of an envelope (see, e.g. ICrosas & MentenI 
Il997h . 

CO is produced abundantly under the thermodynamical 
equilibrium conditions in the atmospheres of such AGB stars, 
even at photospheric temperatures, ~ 2000 K, dTsui ill 1964 ) and 
is consequently observed easily via its near-infrared vibrational- 
rotational transitions. Toward the archetypical, high mass-loss 
C-rich object IRC+10216, Sahai & Wannier (1985) observed 
vibration-rotational e mission and a bsorption in the 4.6 /urn fun- 
damental band, while [Keadv et al.l (IT9881) additionally observed 
absorption from the first overtone band at 2.3 yum. 

Given the extremely high temperature and density in the re- 
gion from which the above transitions are arising, one might also 
expect to find pure rotational emission to be observable from it. 
However, pure rotational transitions of CO in excited vibrational 
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Vibrational State v 

Fig. 1. The open squares with the y-scale to the left show the 
decreasing dipole moment from u = to 4. Between d = 4 and 5 
the dipole moment switches polarity from the negative polarity 
being at the C atom for y < 4 to the negative polarity being at the 
O atom for u > 5. The filled circles show the much larger effect 
in the squares of the dipole moments which are proportional to 
the intensities of the lines. The y-scale is on the right hand side. 



states had not be en reported in the astronomical literature un- 
til very recently (iPatel et al.l (l2009ah . see also below), while ro- 
tational lines of molecules chemically related to CO, although 
with considerably lower vibrational energies, iSyib, have been 
found towji-d IRC+10216. These include CS (Turn er 1987al) and 
SiS (lTurnedll987bl) . which have E^ih = 1 272 cm ' jMiiller et al.l 
l2005h and 744 cm '' (iMiiller et al.ll2007h . respectively; see also 
dPatel et al.ll2009bh . Moreover, vibrationally excited SiO maser 
emission has been detected from the envelopes of many hun- 
dreds of oxygen-rich long period variable stars and red super- 
giants, and b oth its v =1 and 2 (iSyib = 1230 and 2447 cm"', 
respectively (ISanz et al ] l2003l) ^ J = 1-0 lines are observed 
even toward three high-mass protos tellar objects embedded in 
molecular clouds ( IZapata et al.ll2009l) . 

It is noteworthy that for all three of these molecules 
transitions of even higher excited vibrational states and rare 
isotopologues ha ve been o bserved, n amely v — 4 for 
SiO (Cernicharoetai] |1993b and S iS (lAgundez et al.lj2008l: 
IPatel et a l. 2009b) , and t; = 3 for CS ( lAgundez et al.ll2007h and 
^'^SiO jGonzalez-Alfonso et al.lll996h. 

In the case of CO, (IScoville & Solomon|[T978l) reported the 
detection of its y = 1 , 7 = 1 - transition tow ard IRC-H10216. 
Howev er, that "detection" was soon rebutted bv lCummins et al.l 
(Il980l) . They identified the observed spectral feature as due to 
C4H, whose frequency is very close indeed to that of the v - \, 
7=1-0 CO line. Moreover, C4H exhibits generally strong lines 
in this source. It took roughly 30 years until a genuine detection 
of pure rotational CO emission lines {v-\,J-2- \ and 
3-2) with the Submillimeter Array was reporte d (of course in 
IRC H- 102 16) very recently bv lPatel et af] (l2009ah . 

We have recorded rotational transitions of vibrationally ex- 
cited CO belonging to four isotopic species, CO, '"'CO, C'^O, 
and '^C''^0 up to y = 3 in excitation and up to 2 THz in fre- 
quency. In order to facilitate the detection of radio lines of vi- 
brationally excited CO, we present here our results for the main 
CO isotoplogue. The full account of our work, which involves a 
combined fit with a variety of rotational and ro-vibrational data, 
will be presented elsewhere (Gendriesch et al., in preparation). 



2. Experimental details 

Transitions with frequencies below 1 THz have been recorded 
with the Cologne Terahertz Spectrometer (CTS) which has been 
described in detail by Winnewisser et al. (1994). It uses broad- 
band tunable, phase-locked backward-wave oscillators (BWOs) 
as powerful sources, and a magnetically tuned, liquid helium 
cooled hot-electron InSb bolometer as detector The Cologne 
Sideband Spectrometer for Terahertz Application (COSSTA) 
was employed for the measurements between 1.75 and 2.0 THz. 
In this case, the terahertz radiation was generated by mixing of a 
fixed-frequency far-infrared laser and a BWO as tunable side- 
band generator An InSb bolometer was used again as detec- 
tor. Further de tails on this spectrometer system are available in 
iGendriesch et al. (2000). 

The pressure of CO in the 3 and 2 m, respectively, long ab- 
sorption cells was in the 1 - 4 Pa range. The measurements were 
carried out in a static mode or in a slow flow through the cell. 
One has to keep in mind that in the terahertz region a flow speed 
of several m/s causes frequency shifts of a few tens of kilohertz. 

The smaller dipole moment in y = 1 compared with the vi- 
brational ground state reduces the intensities by about 40 %; in 
y = 3 the intensity is decreased by almost one order of magni- 
tude, see Fig. [T] The high excitation energy reduces the intensi- 
ties much more, by about 4.5 orders of magnitude for each vibra- 
tional quantum at 300 K. For these reasons, vibrationally excited 
CO could not be observed even under very favorable conditions 
and with long integration times at room temperature for y > 1 . A 
glow discharge effectively increased the vibrational temperature 
such that lines of vibrationally excited CO could be recorded be- 
low 1 THz with the CTS. The available source power at 2 THz 
is considerably lower with COSSTA. Therefore, a modulation of 
the discharge (on/off) at ~ 5 Hz was employed in addition to the 
usual frequency modulation to reduced baseline effects and to 
improve the spectrometer stabiUty. T he measurements were car- 
ried o ut similarly to those on HNC dThorwirth et al.ll2000h and 
DNC ("Brunk en et ani2006l) in their ground and V2 - 1 excited 
vibrational states. 

A high accuracy of a transition frequency in our measure- 
ments depends very much on the flatness of the baseline, the 
absence of other, interfering lines in the vicinity, and on a 
large signal-to-noise ratio (S/N). Relative accuracies of about 
10 have been reached for strong, isolated lines recorded in 
Doppler mode up to 1 THz as well as near 2 THz fairly com- 
monly as demonstrated in the investigations of, e.g., H2CO 
(Briinken et al. 2003) and SO2 in its y2 = and 1 vibra- 
tional states dMiiller & Briinkeiil 120051) . Under favorable con- 
ditions, better accuracies can be reached even in Doppler- 
limited measurements, as shown in the studies of H CN, y2 = 1 
dThorwirth et al.ll2003l) . and SiS (Muller et al."2007). The study 
of rotational transitions of CO by Varberg & Evenson (1992|) 
demonstrates that such high relative accuracies are possible also 
in other laborato ries and at even h i gher f requencies (beyond 
4 THz); see also IWinnewisser et al.1 dl997h for an updated list 
of these frequencies and uncertainties. It is worthwhile mention- 
ing that in sub-Doppler measurements relative accuracies of bet- 
ter than 10~^ have been reached e. g. in measurements of CO 
dWinnewisser et al.lll997h and '^CO dCazzoli et al.ll2004t) . 



3. Analysis and discussion 

Predictions of the CO, u - 1 to 3 rotational spectrum were 
based initially on infrared mea surements as summarized, e.g. in 
ICoxon & Haiigeorgioiil d2004l) . All but the very weak 7=1-0 
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'C'^O (J = 16 — 15) 



v=3 v=2 v=1 v=0 




1790,97 1790,98 1807,76 1807,77 1824,55 1824,56 1841,34 1841,35 

Frequency / MHz 



Fig. 2. The 7 = 16 «— 15 transitions are shown for v - to 3. 
Note that the signal-to-noise ratio is very good even in the v = 3 
state. The lower state energy of the transition is 6800 cm"', see 
Table [T] The lines appear as second derivatives of a Gaussian 
line-shape because of the 2/-modulation. 

transition were recorded for y = 1 below 1 THz, fewer transitions 
were recorded for higher excited states. Two additional transi- 
tions could be recorded for each vibrational state with COSSTA 
near 2 THz. Fig. |2] displays the very good signal-to-noise ra- 
tios obtained for the y = to 3 states. Altogether between 9 
and 4 experimental transition frequencies were obtained in the 
course of the present investigations. They are given in Table [1] 
together with their uncertainties. Also included are the calcu- 
lated frequencies up to 7 = 20 - 19 with uncertainties, the lower 
state energies, and the line strengths S/u^. 

The spectroscopic parameters have been determined in a 
weighted least-squares fit, and are listed in Table [Altoget her with 
those of the ground vibrational state from .Winnewisser et alj 
( Il997h . The parameters B„ and D„ have been determined very 
well. Even the sextic centrifugal distortion term i/„ has a stan- 
dard deviation that is less than 10% of its value for u = 1, but 
it is barely determined for v - 2 and 3. The rotational constant 
Bu decreases with vibrational excitation, as is generally the case. 
The change in D„ is very small, but still determined with sig- 
nificance; the latter is not the case for //„. The rms error of the 
fits, also known as the weighted standard deviation of the fit, are 
0.561 for V = 1 and 2 suggesting the error estimates of the rest 
frequencies to be conservative. The rms error of 0.2 for v = 3 
is not meaningful because 3 parameters have been determined 
from 4 transition frequencies only. 

Table [U provides calculated rest frequencies up to 2.3 THz. 
The small values of the predicted uncertainties suggest that re- 
liable predictions of the rotational spectrum of CO v = 1 and 
probably 2 can be made up to at least 3 THz. Predictions to 
even higher frequencies will be available in the catalog se ction of 
the Cologne Database for Molecular spectroscop}0 (.Miiller et alj 
l200lLl2005h . 

4. Astronomical outlook 

As mentioned above, pure rotational emission from the ii = 
1,7 = 2-1 and 3-2 lines has very recently been detected 
towards the hot, innermost circumstellar e nvelope of the hig h 
mass-loss carbon AGB star 1RC-H10216 dPatel et all [2009ah . 



Potentially, other vibrationally excited rotational lines from CO 
are detectable toward this and similar objects as well. The obser- 
vations just reported, made with the Submillimeter Array, have 
an angular resolution of a few arc seconds and by no means re- 
solve the emitting region, which is estimated to have more than 
10 times smaller dimensions, ~ !'2 of ~ 30 AU for a distance of 
150 pc (Crosas & Mentenl ll997h . Resolving such small regions 
is outside the capabilities of existing (sub)millimeter interferom- 
eters. 

However, the Atacama Large Millimeter Array (ALMA), 
currently under construction on the ~ 5000 meter high Llano de 
Chajnantor in Chile, with its fifty 12 m diameter antennas will 
have a much higher angular resolution and a much larger collect- 
ing area than all existing instruments. Consequently, ALMA will 
have a vastly larger brightness temperature sensitivity, ATb 
[A^(A^ - l)Aeff]"'(v0B)"^(Av)"'^2. Here is the number of an- 
tennas (currently 50 planned and 64 as a goal), Agff the effective 
area of a single antenna, v the observing frequency, 0^ the syn- 
thesized beam size and Av the frequency resolution. 

Using the online ALMA Sensitivity Calculator TooQ on 
ESO's ALMA WWW site we employ inputs equivalent to v = 
340 GHz (comparable to the frequencies of vibrationally ex- 
cited 7 = 3-2 CO line) and Av = 1 MHz (correspond- 
ing to 0.88 km s"' velocity resolution; i.e., appropriate for the 
~ 7 km s"' width of the lines observed toward lRC-i-10216), we 
calculate for 1 .2 mm precipitable water vapor (mediocre weather 
conditions for the ALMA site) an rms flux density sensitivity of 
3.2 mJy for 1 h of integration time. In an extended configura- 
tion, with a maximum baseline length, Bmax, of 10 km, a beam 
of FWHM Ob - c/(vBmax) radians, or 22 milli arcseconds could 
be synthesized with ALMA (at v = 340 GHz; c is the speed of 
light). With this Ob and the flux density sensitivity quoted above, 
we calculate that an rms brightness temperature sensitivity of 
70 K could be achieved with ALMA at that resolution, which 
corresponds to a quarter of the diameter of lRC-i-10216's radio 
photosphere (Menten et al. in prep.). 

Given that the brightness temperature of the gas emitting vi- 
brationally excited CO is in high likelihood hotter than 1000 K 
(Patel et al. 2009a), spatially resolved imaging will allow studies 
of the distribution and dynamics of the hot material close to the 
stellar atmosphere of lRC+10216 and similar evolved stars. This 
will provide important insight into the origins of the mass loss 
process in such objects. These high resolution observations that 
will also easily image the radio photosp heres of such near-by red 
giant stars (see, e.g., iReid & Mentenll2007l) will be eminently 
feasible since their stellar continuum emission will provide am- 
ple flux for self calibration. 
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Table 1. Measured transition frequencies" (MHz) of vibrationally excited CO, transition frequencies" (MHz) calculated from the 
final set of spectroscopic parameters, lower state energies Eio (cm"') and line strengths Sfj.^ (10""' D^) 
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4603.389 


46.526 


14 




1 695 120.497 (32) 


4656.172 


49.273 


15 


1 807 764.748 (15) 


1807 764.752 (18) 


4712.716 


51.907 


16 


1 920 338.703 (80) 


1920 338.633 (48) 


4773.016 


54.408 


17 




2032 837.75 (11) 


4837.072 


56.797 


18 




2 145 257.73 (21) 


4904.880 


59.032 


19 




2 257 594.17 (33) 


4976.438 


61.146 


!; = 3 







112123.087 (12) 


6350.439 


1.207 


1 




224 241.770 (22) 


6354.179 


2.410 


2 




336 351.644(30) 


6361.659 


3.603 


3 




448 448.307 (35) 


6372.879 


4.781 


4 




560 527.354 (37) 


6387.837 


5.939 


5 


672 584.390(40) 


672 584.384 (37) 


6406.534 


7.075 


6 




784 614.995 (40) 


6428.969 


8.182 


7 


896 614.750(100) 


896 614.785 (49) 


6455.141 


9.255 


8 




1 008 579.357 (65) 


6485.049 


10.290 


9 




1 120 504.311 (85) 


6518.692 


11.284 


10 




1232 385.25 (11) 


6556.068 


12.232 


11 




1344 217.78 (12) 


6597.176 


13.132 


12 




1455 997.51 (13) 


6642.014 


13.976 


13 




1567 720.05 (12) 


6690.581 


14.762 


14 




1679 381.012(94) 


6742.874 


15.488 


15 


1790 976.008 (50) 


1 790 976.004 (49) 


6798.892 


16.151 


16 


1902 500.636(100) 


1902 500.647 (96) 


6858.633 


16.747 


17 




2013 950.56 (25) 


6922.094 


17.275 


18 




2125 321.36(47) 


6989.272 


17.728 


19 




2 236 608.68 (78) 


7060.165 


18.110 



" Numbers in parentheses are one standard deviation in units of the least significant figures. 

Table 2. Spectroscopic parameters" (MHz) of CO in excited vibrational states compared with that in the ground vibrational state 



Parameter 


V = O* 


V = 1 


V = 2 


V = 3 




57 635.968 019 (28) 


57 111.243 12(66) 


56 586.555 00(157) 


56 061.910 5 (59) 


£>„ X 10' 


183.504 89 (16) 


183.4903 (71) 


183.4740(234) 


183.511 (63) 


H„ X 10' 


171.68(10) 


177.9(145) 


163. (53) 


241. (134) 



" Numbers in parentheses a re one standard deviation in units of the least significant figures. 
'' IWinnewisser et al.1 ( Il997i) 



